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SUMMARY 


Existing  work  on  the  pressure  fluctuations  in  turbulent  shear  flows  is 
briefly  reviewed  with  special  reference  tc  the  problem  of  wall  turbulence. 

An  approximate  theory  for  the  pressure  fluctuations  on  the  wall  under 
both  a  turbulent  boundary  layer  and  a  wall  jet  is  given  and  indicates  in 
the  latter  case  an  intensity  many  tir.es  that  corresponding  to  the  flow 
over  a  flat  plate  at  zero  pressure  gradient,  an  typified  by  measurements 
on  the  wall  of  a  wind  tunnel.  Experiments  on  a  wall  jet  confirm  these 
predictions  and  details  of  the  few  preliminary  data  are  presented. 

The  results  from  the  wall  jet  suggest  that  the  intensity  of  the  pres¬ 
sure  fluctuations  in  the  regions  of  adverse  pressure  gradient,  on  wings 
ar.d  bodies,  approaching  and  beyond  separation  will  be  higher  than  in 
regions  of  zero  pressure  gradient. 

Appendices  are  included  which  deal  with  the  necessary  extensions  to 
the  analysis  tc  fit  the  velocity  correlation  functions  as  measured  by 
Grant  (1958),  the  effects  of  time  delay  and  eddy  convection. 


SOMMAIRt 

Les  travaux  effectues  jusqu’ ici  sur  les  variations  de  prussicn 
un  ecculerrer.t  turbulent  er.  cisaj  llerer.t  sent  sermairement  exposes  av.;c 
mention  parti c'iliere  du  problems’  dn  la  turbulence  a  lu  paroi. 

Les  auteurs  presnntor.t  une  th6orio  approch.ee  pour  les  variations  du 
pression  a  la  paroi  sous  V  influence  tant  d’ une  ooucb.e  limite  turbulentu 
£ue  dj  un  jet  d  la  paroi,  ri‘ aprhs  l.vpuelle,  dans  ce  dernier  eus,  les 
variations  du  pression  sor.t,  sensiblement  plus  fortes  cue  cellos  qui 
correspondent  a  1’ ecoulener.t  sur  une  plaque  plane  an  1’ absence  da  gradi¬ 
ents  de  pression,  air.si  que  le  dementrent  les  mesures  eifuctuees  sur  la 
pare:  d' nr.  snufflerie.  Cette  prediction  cst  confirmee  par  les  resultats 
d’ essais  sur  un  jet  a  la  paroi;  les  quelques  resultats  prt?l  imir.aire.s 
obte.nus  sont  preser.tes  ici. 

Ces  resultats  font  per.ser  que  1’  inter.sitc  des  variations  de  pression 
dar.s  les  regions  a  gradients  defavorables,  par  exemple,  sur  les  ailes 
et  les  corps,  quand  1’  ecouler.ent  s’ apprcche  du  point  de  deeclleme-r.t  st 

p?.SS0  ?. U” d?  1 ° ,  ><??***  plus  fr.rt'P  rinp  rl&rs  1  ps  nr.s  R  grarlipri^s  niils, 

Le  rapport  comprend  des  annexes  traitant  des  extensions  a  apportcr  a 
V  analyse  pour  .'a  rendre  applicable  aux  fcr.cticr.s  de  la  correlation  des 
vitesses,  ainsi  qu' elles  or.t  ete  mesurees  par  Grant  en  1958,  a  1' influence 
du  retard  temporel  et  a  la  convection  des  tonrbillcr.s. 
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ON  SURFACE  PRESSURE  FLUCTUATIONS  IN  TURBULENT  BOUNDARY  LAYERS 
G.M,  Lilley  and  T.H.  Hodgson* 


1.  INTRODUCTION 

A  knowledge  of  the  intensity  and  spectra  of  the  pressure  fluctuations  in  turbulent 
shear  flows  is  required  in  a  wide  range  of  aeronautical  and  hydrodynamic  problems 
today.  Such  problems  range  from  aerodynamic  noise  generated  by  turbulent  motion,  the 
vibrations  of  the  aircraft  skin  at  high  speeds  and  the  transmission  of  noise  to  the 
cabin  and  cockpit,  to  mention  just  a  few, 

Although  work  in  this  area  was  pioneered  by  Heisenberg1,  Obukhoff2,  Batchelor3  and 
Kraichnan*’5,  it  is  only  relatively  recently  that  experimental  data  have  been  obtained 
to  check  the  theoretical  data,  and  to  set  the  pattern  for  investigations  into  more 
complicated  situations,  where  the  theory  at  best  would  be  very  tentative.  The  work 
of  Kraichnan  is  of  particular  interest  for  it  deals  with  pressure  fluctuations  in  the 
presence  of  a  mean  shear,  and  also  in  the  case  of  wall  turbulence,  and  therefore  has 
direct  application  to  the  problem  of  wall  pressure  fluctuations  under  a  turbulent 
boundary  layer.  (The  earlier  work  of  Heisenberg,  Obukhoff  and  Batchelor  cpqpidered 
only  the  case  of  isotropic  turbulence).  Kraichnan  showed  that  on  a  wall  /pv(>4  P  u£) 

~  /3cj,  where  (2  is  a  factor  between  2  and  12.  Experimental  results  obtained  by 
Wlllmarth'1  and  Harrison7  confirmed  Kraichnan’ s  predictions  and  gave  values  of  p 
between  2.  5  and  5.  0. 

The  effect  of  Mach  number  on  wall  pressure  fluctuations,  although  of  obvious  current 
importance,  will  not  be  considered  here.  Indeed,  the  flow  will  be  assumed  incompress¬ 
ible  throughout  and  the  problem  of  boundary  layer  noise,  that  is  the  noise  radiated 
away  from  the  surface,  will  hardly  be  touched  upon.  Our  attention  will  mainly  be 
restricted  to  problems  of  wall  turbulence,  including  the  wall  jet,  and  will  not  con¬ 
sider  in  the  same  detail  pressure  fluctuations  in  free  turbulence,  A  review  paper 
covering  the  items  omitted  by  the  authors  would  naturally  be  of  considerable  topical 
interest,  but  it  was  felt  that,  bearing  in  mind  the  considerable  efforts  present  in 
this  area  today  and  the  present  state  of  flux  of  knowledge  in  the  subject,  a  greater 
need  was  for  a  fundamental  appreciation  of  relatively  simple  flow  models.  The  surface 
over  which  the  fluid  flows  will  be  treated  as  rigid,  and  r,o  account  will  be  .tivsn  of 
the  response  of  the  structure  to  pressure  fluctuations. 

The  theory  of  the  pressure  fluctuations  in  wall  turbulence,  including  the  wall  jet, 
will  be  treated  on  similar  lines  to  the  method  used  by  Kraichnan5.  In  this  method  the 
intensity  of  the  pressure  fluctuations  can  be  obtained  once  the  two-point  velocity 
correlations,  mean  velocity  gradient,  and  the  turbulence  intensity  and  scale  are  known. 
Experimental  results  obtained  by  Townsend8,  Laufer5,  and  Grant10  will  be  used  to  find 
p*  on  the  wall.  For  the  wall  jet  the  mean  flow  theory  due  to  Glauert11  will  be 
adopted,  together  with  results  in_the  mean  flow  by  Bakke1'  and  Bradshaw  and  Love13. 

The  prediction  in  this  case  for  p'  will  be  compared  with  those  obtained  from 
measurements. 


'Department  of  Aerodynamics,  College  of  Aeronautics,  Crar.fieid,  England 
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Finally,  in  order  to  avoid  confusion  in  the  theoretical  treatment  and  the  discussion 
of  results,  the  authors  stress  that  they  are  dealing  with  the  problem  of  pressure 
fluctuations  in  a  pseudo-incompressible  turbulent  shear  flow.  They  prefer,  and  agree 
that  it  is  one  of  personal  preference  and  not  one  agreed  by  convention,  to  refer  to 
boundary  layer  noise  as  the  sound  energy  radiated  away  from  the  turbulent  flow.  This 
latter  problem,  considered  by  Curie1",  Phillips15,  and  Doak16,  though  connected,  has 
obvious  differences  from  the  present  treatment. 


2.  THEORY  OF  THE  PRESSURE  FLUCTUATIONS  IN  A  SHEAR  FLOW 


In  most  instances  where  the  pressure  fluctuations  are  significant  the  fluctuations 
in  the  fluid  density  are  significant  also.  The  present  problem  is  no  exception. 
However,  since  we  will  be  more  concerned  with  the  pressure  fluctuations  Inside  the 
turbulent  shear  flow,  which  in  our  problem  will  be  an  essentially  ‘low  speed  flow'  , 
than  with  the  noise  radiated  from  it,  we  may  safely  assume  that  the  flow  is  incompress¬ 
ible.  The  equations  of  continuity  and  motion  are  therefore  respectively 


(1) 


CPU  i 
3 1  3  x 


■  pu,u 


j 


i  j 


3p 

3  X; 


If  we  take  the  divergence  of  (2)  and  use  (1)  we  obtain 
pressure: 


yVUj  (2) 

the  following  equation  for  the 


VJp 


3  Vu< 

p - I— »■ 

Bx^Xi 


(3) 


This  shows  that,  whereas  in  inviscid  steady  flow  the  pressure  at  a  point  follows 
immediately  from  the  dynamic  pressure  at  the  same  point,  the  pressure  at  a  point  in  a 
turbulent  flow,  since  it  obeys  an  equation  of  the  Poisson  type,  is  governed  by  fluctua¬ 
tions  in  velocity  throughout  the  entire  flow  and  not  just  at  the  field  point.  Since 
the  wall  jet  comprises  an  inner  boundary  layer  and  an  outer  jet  mixing  region  we  will,  ' 
for  convenience  in  finding  the  solution  to  (3),  consider  this  more  general  flow  only. 

The  solution  for  the  boundary  layer  on  a  fiat  plate  is  then  a  particular  case. 


In  the  wall  jet  the  mean  pressure  is  approximately  constant  everywhere,  and  the  mean 
velocity  field  is  radially  symmetric  so  that  in  terms  of  cylindrical  polar  co-ordinates 
(R,  0,  Z),  the  velocity  components  are  given  by  uR(R,  Z),  u^  =  0  and  uz(R,  Z) . 

Since  Equation  (3)  above  can  be  written  alternatively  in  vector  notation  as 

V2p  -  -/V-(u-V)u  (3a) 

we  find  cn  expanding  the  right  hand  side  and  making  the  usual  boundary  layer  assumption 
that  (see  Appendix  A) 


3 


(■>  f  "N  —  f  —  \ 

cu:  cUp  duz  oup  ) 
-nr*  -  2  — ^  v^- J 
,3R  or  3r  l'Z j 


where  (u^,  u^,,  u^)  denote  the  fluctuating  components  of  velocity  and  (uR,  u^,  uz)  the 

mean  values. 

3  Up  3uH 

But  — 1  «  — —  so  that  finally  we  have 

3r  3z 

V‘d  =  -  ^P~  (5) 

3r  3z 


where  tr7(R,z) 


is  the  local  nean  shear. 


This  result  is  similar  to  that  found  by  Kraichnan5  for  the  boundary  layer  flow  over 
a  flat  plate  with  a  mean  velocity  in  the  (x,)  direction  varying  with  (x2)  only,  the 
direction  normal  to  the  plate. 

In  order  to  simplify  the  notation,  in  what  follows,  tve  will  drop  the  subscripts  on 
r  and  the  primes  on  uz.  Thus  p  and  u  will  now  refer  to  the  fluctuating  pressure  and 
velocity  respectively,  and  r  is  the  mean  shear.  Therefore  we  can  write  Equation  (5a) 
in  the  equivalent  fore 

V<p  =  -  2pT-— =  =  -  q(x,  t)  (6) 

oR 

where  q(x,  t)  is  the  source  density.  le  are  here  assuming  that  the  pressure  fluctua¬ 
tions  due  to  the  amplifying  effect  of  the  mean  shear  on  the  turbulence  are  greater 
than  that  due  to  the  interaction  of  the  turbulence  on  itself.  Kraichnan5  justifies 
the  use  of  this  assumption  by  showing  that,  approximately,  the  former  will  give  rise 
to  a  root  nean  square  value  of  pressure  about  10  db  higher  over  that  due  to  the  latter. 

The  solution  of  (6)  can  be  found  by  finding  an  appropriate  Green  function,  G(x,  x'), 
which  satisfies  the  boundary  conditions.  Since  the  Green  function  satisfies 

V2o(x,  x')  =  -  5(x  -  x')  (7) 

where  S(y)  is  the  Dirac  delta  function,  the  appropriate  solution  of  (6)  is 


P(x,  t) 


/  q(y,  t)  G(x,  y)  dy  +  /  G—  (y,  t)  -  p(y,  t)  —  dS(y)  (8) 

J  ^  n  *-  A  ri 


where  V  is  the  total  volume  occupied  by  the  flow  and  S  is  the  total  area  of  the  plane 
over  which  the  fluid  flows,  provided  that  at  some  initial  tine  p  and  3p/3t  vanish. 
Here  n  is  the  normal  to  s,  measured  tc  the  surface  from  the  volume. 
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If  we  write  G_  the  Green  function  which  vanishes  on  the  plane  S,  G+  the  Green 
function  given  by 


on  the  plane  S,  and  G,  the  Green  function  for  unbounded  space,  we  have 


(9) 


where  G1  is  the  Green  function  for  the  inage  point  in  the  boundary,  Hence  we  can  find 
three  equi valent  solutions  of  {61,  which  depend  on  the  particular  choice  of  Green 
function.  Thus 


But  on  S  we  have 


on  3n 


and  so  (10)  can  be  written  alternatively 


P(x,  t) 


+  GjldV  * 


Gj)dV  - 


(10a) 


(10b) 


(10c) 


(11) 


(12a) 


(12b) 


(12c) 
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The  particular  choice  of  one  of  these  equations  must  depend  on  what  is  known  about 
p  and  d  p/3  r.  on  the  plane  S.  Clearly  if  cp/cn  -  0  on  S  we  would  choose  (12a). 


Now  the  equation  of  motion  in  the  Z-direction  at  Z  =  0,  where  uR  =  -  uz  =  0.  is 


0 


d  p  c  2uz 
c  z  "  oZJ 


for  the  fluctuating  quantities  and 


0 


(13) 


(14) 


for  the  mean  flow. 


But  from  the  equation  of  continuity 


and  therefore  near  the  wall 


(15) 


4M2\cZ/ 


(IS) 


From  the  measured  distribution  of  u|  near  the  wall  of  a  flat  plate,  channel  or  pipe 
Townsend5  finds  that  (in  our  notation), 


/opY 

\bzj 


_  3 

*  34  *  10"  775 


(17) 


where  rw(R)  is  the  wall  mean  shear  stress  given  by 


(18) 


Although  (cp/oZ) as  shown  by  (17),  is  a  very  small  quantity  we  are  not  immediately 
justified  in  putting  (cp/3z^._Q  =  0  in  (12a). 


However,  if  l,  is  a  typical  length  over  which  the  velocity  is  correlated  and  !„  a 
similar  length  with  respect  to  pressure,  we  see  that  the  contributions  to  p2,  from  the 
volume  ar.d  surface  integrals  in  (12a),  are  in  the  ratio 


100  Z? 


1  2 
l2 


New,  if  anything,  i2  is  less  than  l, 
and  so  we  see  that  the  surface  integral 
Kraichcan' s  approximation5-  Thus  we  r.c 
ar.d  (12c)  are  cf  less  value  than  the  fi 


and  is  certainly  not  large  compared  with 
:r.  (12a)  can  be  neglected  in  agreement  with 
e  that  for  this  problem  the  solutions  (12b) 
•st. 
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We  therefore  write  approximately 

P(x,  t) 


-tL 

2n  J 


x  -  yl  1  -  :x  -  y*|  *)  T(y)  —±.  (y)  dy 

3yR 


(19) 


on  inserting  q("x,  t)  from  (6)  into  (12a),  and  noting  that 

’  /I*  -  _  j_  N', 

4t -  \lx  -  y|  |x  -  y* 1  J 


Gi  =  “ 


(20) 


where  y  =  (yR,  y^,  yz);  y*  s  (yR,  y^,  -  yz),  and  the  volume  of  integration  is  over 
the  half-space  Z  >  0. 

On  the  surface  of  the  plane  x  =  (R,  <p,  0)  so  that 


P(x,  t)z=D 


(21) 


and  the  pressure  covariance,  for  zero  time  delay  between  the  pressures  on  the  plane  at 
x  and  respectively,  is 


P(X,  X  )£„(.  _2 


dy  dz 


(22) 


where  the  velocity  covariance  is 


Rj  i  =  Uj  ■  (y)  u;  (z) 

and  Z  -  (ZR,  zz>- 

If  the  turbulence  is  spatially  locally  homogeneous  we  can  put 


R77  (y  r  Z)  _  — , 

; _  -I - f _  J  -  R77  (r) 

slH  ^  V  uz  (z) 

where  r  -  (z  -  y)  and  R^  is  the  velocity  correlation  coefficient. 


(23) 


Now  the  integral  in  Equation  (22)  with  (23)  is  difficult  to  evaluate  for  general 
values  of  r  and  uz.  However,  the  contributions  to  the  integral  beyond  a  certain  dis¬ 
tance  of  the  order  of  the  length  scale  of  the  turbulence  will  be  sufficiently  small  so 
that  r  and  uz  car.  be  regarded  as  functions  of  jx|  and  yz  only.  This  assumption,  means 

that  so  far  as  the  integral  in  (22)  is  concerned  r  and  uz  vary  only  with  distance 
normal  to  the  plane  and,  in  view  of  the  locally  homogeneous  assumption  in  (23),  the 
integral  can  be  evaluated  as  if  the  flow  were  parallel.  For  convenience  in  its 
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evaluation  we  will  use  a  rectangular  cartesian  tensor  notation  with  yl  and  y3  taken  in 
the  plane  of  the  wail  and  y2  normal  to  the  wall.  Thus  (22)  becomes 


P(?;  x')z«c 


X 


-cc 


X 


f 


03 


X 


X 


yj 

-» 


r 


x 

f  "(x;y2)  r(x;  y2  *  r2) 
|  drs  x  -  y|  |x'  -  y  -  ? 


ui!  (x:y2)  J  u2  (x;y2  +  rj 


R,, 


-  (r)  (24) 


o  r‘ 


If  R22  and  its  derivatives  with  respect  to  vanish  at  infinity  an  integration 
by  parts  of  (24)  leads  to 


P(x;  x')z=c 


(r) 


T(x;  y,)  -r(x;  y  ;  y,  '  r  ) 


(25) 


Now  it  rust  be  noted  that  in  this  analysis  stationary,  and  not  roving  co-ordinates 
are  being  used.  Also,  as  stated  above,  P(x;  x')(  the  pressure  covariance  on  the  plane 
Z  -  0,  is  that  for  zero  time  delay  and  involves  the  instantaneous  product  of  the 
fluctuating  pressures  at  the  points  x  and  x'  respectively.  However,  if  the  frequency 
spectrum  of  the  wall  pressure  fluctuations  is  required  it  is  necessary  to  consider 
both  the  tire  delay  and  the  turbulence  in  a  roving  frame  of  reference  (see  Appendices 
D  and  E). 


Equation  (25)  is  the  important  equation  in  this  paper.  It  shows  that  the  pressure 
covariance  at  a  point  on  the  wall,  due  to  shear  flow  turbulence  docinated  by  the  mean 
flow  shear  (r) ,  is  governed  by  R  the  two-point  ‘lateral’  velocity  correlation 
coefficient,  the  mean  shear  and  the  turbulent  intensity. 


In  Equation  (25)  the  integration  over  y  can  only  be  effected  when  the  values  of 
t/oI  are  known.  Now  we  car.  easily  see  that  the  greatest  values  of  r/ui  exist  in  two 

c  i 

quite  distinct  regions  of  the  wall  jet  at  any  value  of  the  radius  R,  These  regions 
are  the  constant  stress  region  of  the  boundary  layer  like  flow  and  the  middle  of  the 
outer  mixing  region  respectively.  On  the  assumption  that  the  turbulent  inner  region 
of  the  wall  jet  possesses  the  same  structural  similarity  as  the  irxer  region  of  a  flat 
plate,  channel  or  pipe  flow  we  have  that  the  mean  velocity  distribution  in  this  region 
(known  as  the  'Law  of  the  Wall’  )  is 


.  f i ZuJ\ 

In  !  — j 

V  v  J 


A 


K 


(26) 


a 


Bu 

and  r  —  —  =  -21  (27) 

3  Z  KZ 

where  uT  =  Aw/p  is  the  shear  velocity. 

Equation  (26)  is  known  to  hold  for  ZuT/V  >  30  and  Z/Z^  <  0.2,  where  is  the 

thickness  of  the  shear  layer  up  to  the  maximum  velocity  u  =  u  .  K,  known  as  the 
von  Karman  constant,  has  a  value  roughly  equal  to  0.4.  Since  the  motion  in  this  inner 
region  must  be  described  by  universal  functions  of  the  wall  shear  stress,  r  and  the 
viscosity,  the  turbulent  intensity  u^  must  be  proportional  to  u2.  Prom  the  measure¬ 
ments  of  Laufer*  in  a  circular  pipe  it  is  found  that 


(28) 


a  constant  over  most  of  the  inner  region,  (a  is  of  the  order  of  unity).  From 
Equations  (27)  and  (28)  we  see  that 

/T  a  ut 

t/u|  = - — —  (29) 

z  K>  (ZuT/V) 

for  values  of  (Zuy/v)  >  30,  Nearer  the  wall  it  is  found  that 


(30) 


and  it  is  interesting  to  note  that  the  maximum  value  of  r/ajf  occurs  Just  outside  the 
laminar  sub-layer.  In  fact  it  is  very  close  to  the  region  where  the  total  turbulent 
intensity  (u^  -  -  u|)  reaches  its  maximum  value  (Fig.l). 


It  is  found  from  (29),  (30)  that 


f  dZ* 

f(Z*)  3  -  „  •  0.43  (31) 

where  f(Z*)  =  r2  u|  and  Z*  =  (Z  v>r/v).  Thus  the  average  value  of  r 2  u|  in  the 

region  0  <  Z  <  Z^  is  0,43  u*/v^. 

On  the  other  hand  it  is  shown  in  Section  8  that  the  value  of  the  mean  shear  in  the 
outer  mixing  region  is 


-  0.616  u._ 
-  Zm> 


(32) 


where  Zw,  z_  are  the  distances  to  the  points  of  half  and  maximum  velocity  respectively. 
On  the  assumption  that  in  this  latter  region  /u|  -  0.15  u  we  find  that 
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0.0085  u* 


(33) 


This  core  or  less  completes  the  formal  treatment  of  the  pressure  fluctuations  in  a 
wall  jet.  Since  the  contributions  to  p7  from  the  inner  region  and  the  outer  mixing 
region  are  very  different  in  magnitude  (see  (31)  and  (33))  we  will  consider  them 
separately. 


3.  PEES  SIRE  FLUCTUATIONS  DUE  TO  THE  BOUNDARY  LAYER 
TYPE  SHEAR  FLOW 

Proa  (25)  we  fnund  that 


p(x;  x')z.0  = 


oc  a 

D  ( 

f 

33  CD  C 

r  r 

D  CC 

r  r 

4y . 

I  i72 

1  dy3  dr.  j 

dr2 

CL  ' 

\ 

3  -< 

L  -i  A 

J  „  -OD 

2 

dr,  R  (r) 


t(x;  y2)  r(x;  y.  -  r2)J  u2  (x;  y.)J  u2  (x;  +  r2) 


(34) 


JL  f  i  \ 

By.  Vjx  -  y  !y  B  z .  \?'  -  z\j 


If  the  turbulence  is  assumed  isotropic 

f  rf'\ 


R:  ,  (r)  -  If  +  -  )  i - rs  r, 

J  t  y  2  J  J  ‘  2r  J  1 


(35) 


where  f  is  the  longitudinal  velocity  correlation  coefficient  and  f;  =  df/dr.  Hence 

rf'  f' 

R..  (r)  =  f - - r I  (36) 

22  2  2r  2 

This  form  for  R22  does  not  agree  too  well  with  the  experimental  results  of  Grant10 
and  others.  However,  the  errors  involved  in  the  use  of  (35)  can  be  shown  to  be  small. 

If  we  use  the  result  from  Appendix  B  (Eqn.  B - 1 2>  we  find  that  approximately 


P(x,  o)z_ 


8a 

15 


— pi  <  r'  u?  >  J  rf  dr 


(37) 


where  a  is  an  anisotropy  factor  which  is  assumed  to  have  a  value  of  about  1/3,  and 
from  (31) 


(38) 
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If  we  now  insert 


r  =  *,p(w 


into  (37) 


f  dr  =  L 


P(x.  °)z  =  -  1 

- %  0.28  — 

H  P  ui  Z„ 


Zm 

El  T 


cf  %  3cj 


where  we  have  used’  Z  ■=  0.3  Zc;  Z*  uT/V  =  1600  and  c{  =  tw/%  p  is  the  local  skin 
friction  coefficient.  This  result  is  similar  to  that  found  by  Kraichnan5  for  the 
boundary  layer  on  a  flat  plate,  except  that  our  factor  3  replaces  the  range  2  to  12 
as  given  by  Kraichnan  although  we  note  that  strictly  /P(o)/H  P  u*  is  proportional  to 
Cf'*'  and  v-X  ujv. 


Equation  (41)  can  be  compared  with  the  results  of  various  experimenters  for  ordinary 
boundary  layer  flow.  The  measurenents  of  Harrison7  suggest  a  value  of 


while  Willmarttr7  gives 


yfP/iipu^ 


%  0. 0095  =  4.  8  cf 


0.006  =  2.5  cf 


Our  measurement  on  the  wall  of  a  wind  tunnel  with  a  microphone  0.14  in.  diameter 
and  a  boundary  layer  displacement  thickness  of  0.238  in.  gave 


^P^4pu|  =  0. 


008  =  3,8  c, 


All  these  results  are  qualitatively  in  fair  agreement  with  the  theory  and  indicate 
that  in  a  turbulent  boundary  layer  on  a  flat  plate,  in  zero  pressure  gradient,  the 
wall  pressure  fluctuation  arises  mainly  from  disturbances  in  the  constant  stress 
region,  Since  cf  changes  only  slowly  with  distance  we  see  that  /p7  will  only  change 
slowly  with  increase  in  x. 


Finally  we  should  note  that  our  analysis  leading  to  Equation  (41)  is  only  approx¬ 
imate.  Nevertheless  the  more  accurate  analysis  given  by  Lilley18  leads  to  a  value  of 
3. 1  c  t . 


•This  value  is  obtained  from  the  work  cf  GrantlJ  by  noting  the  zero  point  value  of  R;2  (r,  0,  0) 
in  the  cuter  regies  of  a  boundary  layer 


ii 


4.  PRESSURE  FLUCTUATIONS  DIE  TO  THE  'JET  MIXING  REGION’ 

SHEAR  FLO* 

In  the  outer  region  of  the  wall  jet  the  structure  of  the  turbulence  is  largely 
uninfluenced  by  the  wall.  Nevertheless  the  solution  of  Equation  (6)  is  still  given 
by  (19),  for  the  contribution  from  the  *imag_e’  Green  function  must  still  be  included. 
However,  the  problem  is  simpler  because  r*uz  is  approximately  constant  over  the  region 
centred  about  the  middle  of  the  mixing  region  at  any  radius.  Equation  (25)  therefore 
reduces  to 


P(x,  x’) 


X  X 

r 


Z=C 


'■(X)  u f ( x )  J  dy.  I  dy2  J  dy,  J  dr.  j 

-x  ^3  -X  -  x  -  y 


« 

r 

dr.  I  dr 


R22(r) 


-  (  1  \  B  f  1  ^ 

r.  \!x  -  y'fbT,  '\\x!  -  z\  J 


If  again  put  (xf  -  x)  =  0  then  fro c  Appendix  B 


P { X ,  X ) 2=  q 


a; 

82  .  .  —  _  f 

—  p  t‘(x)  u'(x)  rf 
15  2  . 


dr 


where  a  is  the  anisotropy  factor. 

If,  further,  f  -  exp(-r /lQ) 

then  P(x;  x)z=c  = 

On  inserting  the  values  for  r'  u‘  from  (33)  and  putting* 

=  <Z*  -  Zn> 


(45) 


(46) 


(47) 

(48) 


(49) 


then 


! \  p  ul 


**  0. 1 


(50) 


The  value  obtained  from  our  measurements  is 


7 - t  =  0.11  (51) 

n  p  ui 


and  indicates  that  the  ccr.trihut icn 
times  that  from  the  inner  boundary 


layer  like  f 


the 

1  rjj 

■*.  O  w  . 


outer  mixing  region  is  some  ten 
This  result  is  not  surprising 


•This  value  of  I.  corresponds  to  the  half -width 
this  value  seems  unduly  high  but  it  is  in  keep 
ir.  Appendix  C  ar.d  the  measured  spectra. 


of  the  outer  mixing  region.  At  first  sight 
ng  with  the  results  cf  the  large  eddy  analysis 
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■hen  one  remembers  that  the  small  shear  in  the  outer  mixing  region  is  acting  over  a 
distance  many  times  that  of  the  entire  inner  boundary  layer  thickness. 

The  theoretical  value  of  is  independent  of  the  radius  (R)  and  this  is 

confirmed  by  our  experimental  results.  This  altogether  interesting  result  is  we  feel 
most  important,  for  it  shows  that  even  in  a  region  of  flow  where  the  velocity  changes 
in  the  mainstream  direction  are  large  (uT  ^  1/R),  the  pressure  fluctuations  are  still 
largely  determined  by  local  conditions  only.  This  agrees  with  the  results  of  the 
mean  shear  stress  on  the  wall  by  Bradshaw  and  Love13  (Fig, 2).  They  showed  that  the 
skin  friction  coefficient,  although  slightly  higher  than  for  a  flat  plate,  changed 
very  slowly  with  increase  in  radius  beyond  an  inner  radius  where  the  sail  jet  was 
being  established. 


3.  THE  SPECTRIW  AND  SPACE  PBESSUBE  CORRELATION  FOR 
ZERO  TIME  DELAY 

In  Sections  (3)  and  (4)  the  formal  treatment  of  the  pressure  covariance  with  zero 
time  delay  is  given,  but  only  the  results  for  the  mean  square  of  the  pressure  are 
evaluated  in  full.  However,  Lilley?  has  already  evaluated  the  pressure  covariunce  in 
free  shear  flow  turbulence,  and  as  we  have  seen  above  the  results  for  wall  turbulence 
will  be  similar  apart  from  a  numerical  factor.  From  Lilley’s  results  for  P(x;  x') 
with  f  -  exp(-  c2  r2),  we  find  (Fig. 3)  that  P(i;  £  0,  0)  has  a  large  negative  loop 

with  a  2ero  crossing  point  of  org  %  l.  On  the  other  hand  P(x;  0,  0,  |3)  is,  positive 
for  all  values  of  and  falls  to  zero  as  x  much  core  slowly  than  f (r),  as  r  -  ~. 

This  isotropic  turbulence  model  used  in  evaluating  these  results  is  a  little  too  crude 
to  make  comparison  with  experiment  justified,  apart  from  an  order  of  magnitude  basis, 
yet  it  is  interesting  to  note  that  Harrison7  found  similarly  positive  values  for 
P(x;  0,  0,  ~3)„  although 'the  fall-off  at  large  values  of  ~3  was  slower  than  in  our 
results.  However,  this  is  what  one  might  expect  from  the  anisotropy  in  the  large  scale 
turbulence  modifying  the  form  of  R22(r)  at  large  values  of  ir], 

The  wall  pressure  spectrum  function  has  been  evaluated  in  Appendix  C  taking  into 
account  the  large  eddy  structure,  and  in  Appendix  E  allowing  for  convection.  In  both 
these  cases  and  also  in  the  free  shear  flow  turbulence  example,  the  spectrum  at  low 
frequencies  obeys  the  law  Tie.exp(-  Td )  -  (Fig. 4).  In  our  wall  jet  experiments  we  find 
such  a  rise  at  low  frequencies,  in  contrast  with  all  measurements  made  in  pipes  and  on 
tunnel  wall  boundary  layers,  where  the  flatness  of  the  moderately  low  frequency  end  of 
the  spectrum  is  marked. 

Harrison7  has  also  remarked  on  this  flatness  of  the  spectra  at  low  frequencies  and 
has  suggs, stsd!  t h sit  6iC6pt  at  the  very  lowest  frequencies  outside  the  range  of  the 
measuring  equipment,  this  might  be  explained  in  terms  of  the  intermittency  of  the 
boundary  layer.  However,  this  does  not  in  itself  explain  the  relatively  high  energy 
content,  in  the  lower  frequencies. 

A  complete  explanation  of  this  phenomenon  has  not  yet  been  found,  assuming  of 
course  that  the  flatness  in  the  low  frequency  end  of  the  spectra  is  net  associated  with 
spurious  wind  tunnel  effects,  such  as  fan  noise,  flow  noise  extra  and  above  the  boundary 
layer  noise,  or  tunnel  circuit  resonance.  However,  a  clue  tc  the  explanation  may  come 
from,  an  analysis  of  results  similar  to  the  space-time  correlations  of  wall  pressure 
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obtained  by  Willmarth6*.  These  data  show,  as  explained  in  Appendix  D,  that  the  spatial 
pressure  correlation  for  optimum  time  delay  (i.e.  roughly  the  autocorrelation  in  axes 
moving  with  the  mean  convection  speed  of  the  eddies)  does  not  fall  to  negative  values 
at  large  separation  distances  in  contrast  to  the  autocorrelation  measured  in  axes  fixed 
in  the  wall,  This  means  that  the  life  times  of  the  big  eddies  are  being  extended  on 
account  of  the  growing  scale  associated  with  the  slow  increase  in  the  boundary  layer 
thickness.  This  consequent  modification  of  the  autocorrelation  in  moving  axes  at 
large  times  could  produce  some  lift  to  the  low  frequency  end  of  the  spectra  although  a 
fall  off  like  £2  might  still  be  expected  at  the  very  lowest  frequencies. 

On  the  other  hand  it  would  appear  that  in  the  case  of  the  wall  jet,  since  the 
Increase  in  shear  layer  thickness  is  relatively  large  coupled  with  the  rapid  fall  off 
in  velocity  with  increase  in  radius,  the  life  time  of  the  eddy  could  not  be  extended 
in  this  way.  These  differences  between  the  theoretical  and  experimental  spectra  for 
both  the  wall  jet  and  boundary  layer  are  being  further  investigated. 


6.  THE  WALL  JET  -  THEORY  OF  THE  MEAN  FLOW 

The  wall  Jet  has  been  studied  theoretically  by  Glauert11  and  experimentally  by 
Bakke*2  and  Bradshaw3.  Although  Glauert  found  solutions  for  both  the  laminar  and 
turbulent  problems  only  the  turbulent  case  will  be  required  here.  He  found  that  the 
velocity  distribution  in  the  vicinity  of  the  wall  was  similar  to  that  in  a  boundary 
layer  with  zero  pressure  gradient.  In  the  inner  region  one  would  expect  the  ‘law  of 
the  wall'  to  apply  but  for  simplicity  Glauert  assumed  the  Blasius  distribution 


where  uT  is  the  shear  velocity  given  by  uT  =  /(r  /P). 

In  the  outer  region  Glauert  evaluated  the  velocity  distribution  numerically  and 
showed  that  it  is  slightly  fuller  than  that  given  by  a  (l/7)th  power  law.  The  maximum 
velocity  occurs  at  Z/Z^  =  0.125,  for  a  value  of  (uT/V) (Z^  -  Z^)  -  4  x  10\  where  is 
the  ordinate  to  the  point  in  the  outer  mixing  region  for  which  uR/un  =  0.5.  The  value 
Z/Z^  =  0.125,  where  u  -  uE,  differs  from  that  obtained  experimentally  by  Bakke,  but 
this  can  be  explained  from  Glauert' s  analysis,  for  Bakke’ s  value  of(um/iO(Zli  -  Zn)  = 
3.5  x  1 0 3 ,  compared  with  our  value  above  of  4  x  101*. 

The  distribution  of  maximum  velocity  with  radial  distance  follows  the  law 

1 

UR  ~  (53) 

*  R 

at  the  Reynolds  numbers  of  the  tests  reported  here. 


•Results  analcguous  tc  these  have  been  obtained  by  Favre  and  his  co-workers^c  for  the  two-point 
velocity  correlations  with  separation  in  the  streamwise  direction. 


fi?  H '"f '  •»  y  V<mfrp  "  'p-  '  WWfPTW't'’' i  «*•<»«*•  •  ■»♦  *»••••«  -up* -  -  *v»w  w  -»t**  -• '  t’«"M  »W» W  *T 
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The  velocity  distribution  in  the  outer  mixing  region  is  found  by  Glauert  to  be 
given  approximately  by 


where  ~  R° • 9 . 


-*  =  sech^  0.875  1' 


(z  -  \ 
Vzs  -  V 


(54) 


The  experiments  of  Bakke  confirm  these  predictions.  Bradshaw  and  Love  give  some 
results  for  the  variation  of  wall  shear  stress  under  the  wall  Jet.  These  show  that  in 
the  region  of  fully  developed  flow  the  wall  shear  stress  is  some  25%  higher  than  in 
the  corresponding  case  (i.e.  equal  ur  and  Z^)  for  the  flow  over  a  flat  plate  in  zero 
pressure  gradient. 


7.  APPARATUS 

Wall- Jet 

The  test  rig  shown  in  Figure  5  was  geometrically  similar  to  that  used  by  Bakke*2 
for  confirmation  of  the  theoretical  results  obtained  by  Glauert11.  The  jet  was  1.5 
inches  diameter  and  air  was  supplied  from  pressurised  reservoir  tanks  via  a  6  inch 
diameter  throttling  valve  and  200  feet  of  6  inch  diameter  pipe,  part  of  the  laboratory 
ring-main  supply  which  passes  horizontally  18  feet  above  the  test  site,  The  vertical 
downpipe  was  3  inches  diameter  and  10  feet  long  connected  to  the  6  inch  main  through 
a  3  inch  isolating  valve.  This  was  followed  by  a  smooth  contraction  containing  a  wire 
gauze  and  then  5  feet  of  1.5  inch  diameter  .smooth-bore  pipe  fitted  with  a  6  inch 
diameter  flange  forming  the  jet.  The  jet  flange  was  0.75  inch  above  the  plate. 

Two  plates,  each  measuring  4  feet  square,  were  used.  One  was  constructed 
of  light  acoustic  boarding  2  inches  thick  lined  with  thin  Bakelite  sheet;  the  other 
was  of  2  inch  thick  Tufnol  about  100  lb  in  weight. 

The  test  procedure  was  tc  open  the  3  inch  valve  fully  and  to  control  the  air-flcw 
by  the  main  6  inch  valve  at  the  reservoir  tanks.  This  was  necessary  to  keep  the  valve 
and  pipe  noise  *o  a  minimum. 


15 


The  measuring  region  was  3-10  jet  diameters  from  the  jet  axis  over  which  the  jet 
half  width  (distance  to  half-maximum  velocity)  varied  from  0.4  inch  to  1.2  inches 
with  maximum  velocities  of  140  f/s  to  40  f/s. 

Pressure  Transducers 

Ammonium  di-hydrogen  phosphate  crystal  transducers,  type  H213  and  M141,  were  used. 
These  are  manufactured  by  the  American  Massa  Company. 

The  M213  has  an  outside  diameter  of  0.22  inch  and  a  diaphragm  diameter  of  0.14 
inch.  Its  capacity  is  12  pf  and  has  a  level  frequency  response  to  120  kc/'s,  It  is 
used  with  a  low  noise  and  low  microphone  cathode-follower  connected  through  6  inches  of 
low- noise  cable.  The  input  impedance  of  the  cathode-follower  is  200  megohms  shunted 
by  6  pf.  allowing  a  response  down  to  20  c/s.  The  cathode-follower  noise  level  is  5  mV, 
although  below  300  c/s  the  noise  from  the  transducer  is  25  mV.  A  high-pass  filter 
with  a  cut-off  at  250  c/s  can  be  inserted,  but  is  not  normally  necessary.  The  sensi¬ 
tivity  at  the  cathode-follower  in  the  output  terminals  is  -111  decibels  re  1  volt/ 
microbar  (approximately  2.8  microvolt/(dyne/cm2) ) . 

The  signal  was  amplified  by  a  battery  powered  low-noise  amplifier  of  28  decibels 
voltage  gain,  followed  by  an  amplifier  of  94  decibels  voltage  gain.  Root  mean  square 
readings  mere  measured  on  a  meter  of  the  linear  averaging  type.  The  bandwidth  of  the 
amplifiers  was  5  c/s  to  above  500  kc/s. 

The  Ml 41  has  a  diameter  of  0.6  inch,  a  capacity  of  110  pf  and  a  sensitivity  at  the 
cathode-follower  output  terminals  of  -94  decibels  re  1  volt/microbar  (approximately  20 
microvolts/(dyr.e/cm2)).  The  low  frequency  noise  of  this  transducer  is  much  lower  than 
that  of  the  M213  due  to  its  higher  capacitance,  and  the  usable  frequency  response  is 
20  c/s  -  30  kc/s. 

The  transducers  were  mounted  in  the 'light  acoustic-board  plate  in  soft  rubber 
sleeves  and  in  the  Tufnol  plate  with  a  heavy  brass  body  and  ‘0’  ring  suspension,  similar 
to  that  of  Willmarth"7.  The  use  of  the  two  methods  of  fixing  in  the  two  very  different 
plates  was  to  ensure  that  the  measurements  were  not  affected  by  the  vibration  of  the 
transducers  due  to  the  impingement  of  the  jet  on  the  plate. 

Spectrum  measurements  were  made  using  a  set  of  third  octave  filters,  covering  the 
range  40  c/s  -  20  kc/s. 

Ccrre later 

The  correlator  used  in  the  space  and  time  correlations  of  pressure  and  velocity 
fluctuations  worked  on  the  analogue  principle.  The  multiplier  was  of  the  quarter- 
squaring  type  and  used  two  special  squaring  valves.  The  design  was  based  on  that  of 
Miller,  Scltes  and  Scott21.  However,  much  improved  circuitry  was  employed  so  that,  if 
necessary,  multiplication  could  be  accomplished  over  the  bandwidth  DC  -  200  kc/s.  The 
accuracy  was  not  of  the  standard  associated  with  computing  multipliers,  but  is  of  the 
order  of  1  -  2f;  which  is  quite  acceptable  for  the  present  purpose.  The  output  was 
read  on  a  DC  galvanometer  with  a  time  constant  variable  between  2-10  seconds. 
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A  time  delay  is  accomplished  on  a  twin  channel  tape  recorder  of  extremely  low  wow 
and  flutter,  using  one  fixed  and  one  novice  head.  The  recorder  used  0.5  inch  wide  tape 
running  at  75  inch/second. 

The  two  signals  to  be  correlated  were  firet  recorded  on  a  loop  of  tape  giving  a 
sample  length  of  2  seconds  (this  sample  length  could  be  Increased  if  extra  idling 
pulleys  were  fitted  to  the  machine)  and  then  played  back  inserting  the  required  time 
delay  by  manual  control  of  the  moving  head.  The  signals  could  either  be1  recorded  on 
a  frequency  modulation  system  bandwidth  DC  -  20  fee/s  or  on  an  amplitude  modulation 
system  with  a  bandwidth  of  250  c/s  -  100  kc/s.  The  signal  to  noise  ratio  was  42  decibels. 

The  correlator  was  capable  of  correlating  a  100  kc/s  sine-wave  with  15  points  to  the 
period.  The  maximum  time  delay  possible  was  40  millisecs  and  the  minimum  increment  was 
0.  67  .microseconds. 


8.  PRELIMINARY  RESULTS  AND  DISCUSSION 

The  measurements  of  mean  velocity  across  the  wall  Jet  from  3  to  10  diameters  from 
the  Jet  axis  showed  good  agreement  with  the  similarity  profiles  obtained  theoretically 
by  Glauert11  and  experimentally  by  Bakke12  (Figs. 6  and  7).  The  maximum  velocity  varied 
as  R"*’55,  where  R  is  the  distance  from  the  Jet  axis,  and  the  jet  half  thickness  Z« 

A  fT  ."X 

varied  as  R  *  ,  The  maximum,  velocity  occurred  at  Z/Z^  =  0.15  corresponding  to  a 
Reynolds  number  uc(Z^  -  Zn)/v  =  2  x  10\ 

The  measurements  of  the  wall  pressure  fluctuations  from  3  to  10  diameters  from  the 
jet  axis  showed  that 

Jt?  kpv-l  -  0.11 

and  was  approximately  independent  of  radius  (Fig, 8).  The  corresponding  spectra 
(Fig. 9)  obeyed  a  similarity  law  on  the  basis  of  the  frequency  parameter  co(2^  -  Z^. ) /u^, 

A  typical  autocorrelation  is  shown  in  Figure  10. 

Comparative  results  were  also  obtained  in  the  College  of  Aeronautics  20  in,  x  li  in. 
low-turbulence  wind  tunnel  at  a  position  in  the  working  section  where  the  displacement 
thickness  was  0.286  in,  at  a  freestreaa  speed  of  126  f/s.  The  corresponding  value  of 
the  skin  friction  coefficient,  c,,  obtained  from  the  measured  velocity  profile,  was 
0.  0022-  The  value  of  y/pT/fep u*  was  0,008  which  is  in  fair  agreement  with  the  measure¬ 
ments  of  Willmarth  and  Harrison.  However,  this  result  cannot  he  relied  upon  quantita¬ 
tively  because  in  this  tunnel  the  low  frequency  extraneous  noise  level  is  not  low, 
although  inevitably  its  contribution  to  the  total  pressure  energy  is  small.  Perhaps 
what  is  more  important  is  that  it  provides  a  check  both  on  the  accuracy  of  the  instru¬ 
mentation  used,  and  or.  the  values  of/pT/^pu^  as  found  in  the  wall  jet  using  the  sane 
pressure  transducers. 

The  large  values  of  /p^spu^  as  measured  on.  the  wall  jet  are  not  special  to  that 
case  and  have  been  found  by  Qwen“  to  exist  on  wings  in  regions  of  separated  flow, 
However,  the  experimental  data  on  the  effect  of  pressure  gradient  on  wing  pressure 
fluctuations  are  sc  few  that  we  cm  only  speculate  that,  on  the  basis  of  our  theory, 
the  surface  pressure  fluctuations  will  be  high  when  either  the  value  of  c*.  Is  high, 
say  at  transition,  or  the  mean  shear  layer  is  very  thick,  say  approaching  and  beyond 
separation. 
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9.  CONCLUSIONS 

A  review  of  theoretical  and  experimental  work  on  wall  pressure  fluctuations  in 
turbulent  boundary  layers  has  been  presented.  The  theory  due  to  Kraichnan  is  modified 
and  extended  to  include  the  separate  effects  of  the  large  eddy  structure  and  the  con¬ 
vection  of  the  eddies,  and  the  treatment  covers  both  the  turbulent  boundary  layer  on 
a  flat  plate  and  the  wall  jet.  The  latter  case  is  presented  for  it  provides  data  in 
the  important  practical  case  when  the  turbulence  is  being  subjected  to  a  rapid 
variation  in  mean  shear,  not  unlike  that  associated  with  the  flow  approaching  and 
beyond  separation, 

Preliminary  theoretical  and  experimental  data  for  the  wall  jet  give  values  of 

JiFPbP  u*  many  times  that  of  corresponding  measurements  on  a  flat  plate  with  zero 
pressure  gradient.  These  results  are  explained  in  terms  of  the  greater  thickness  of 
the  shear  layer,  the  relatively  high  intensity  of  the  turbulence,  and  the  presence  of 
relatively  large  eddies  in  the  flow. 

The  present  work  is  intended  as  a  basiB  for  further  work  at  both  low  speeds  and  in 
the  more  important  practical  areas  of  supersonic  and  hypersonic  flight. 
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SMALL  CDOY  09  (l-y^)  «*p  (-  <//)  4/ 0,10 


Fig. 11  Comparison  between  theory  and  experiment  for  RJ?(r,,0,0)  in  a  turbulent 

boundary  layer  y/S 0  =  0.076 


LASGf  EDDY  O'!  txp  *,  * 

SMALL  EDDY  09  txp  (-rJ/j) 


Fig. 12  Comparison  between  theory  and  experiment  for  R22(0,  r 2,  0)  in  a  turbulent 

boundary  layer  y/50  =  0.  053 
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SMALL  EDDY  CW  (l- tip  (- Y^) 

Pig. 13  Comparison  between  theory  and  experiment  for  R^fO, 0, r3)  in  a  turbulent 

boundary  layer  y/&c  =  0.052 


Comparison  between  autocorrelations  of  R.,(r  ,0, 0;  t)  in  the  fixed  and  moving 

frames  in  grid  turbulence 


Fig. 14 
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P  (rr\fcT,r,o,o/r)  =(l“<r(r1-Vfc'n)(l-V|T|)  exp  (-<r  |  r,- VCT  l)  exp(-V  ITI) 
WITH  v-<rvl  vc  =  0-0  V*  O' 033  ft  V'/vc{5sI/jo 

Fig.  15  Boundary  layer  wall  pressure  fluctuations  space-time  correlation  - 
Comparison  between  theory  and  experiment 


APPENDIX  A 


The  Pressure  Equation  for  the  Vail  Jet 


The  equations  of  continuity  and  notion  for  an  incompressible  flow  are,  in  vector 
notation: 


7-  (u)  =  0 


(A-l) 
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where 
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If  in  the  cylindrical  co-ordinate  system  (R,  0,  Z)  we  write  the  corresponding 
velocity  components  (uR,  u^,  uz),  then  the  pressure  distribution  equation,  found  from 
subtracting  the  time  derivative  of  (A-l )  from  the  divergence  of  (A-2), 
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can  be  written  in  full  as 
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Now  the  mean  velocity  u  -  (uR,  0,  uz)  satisfies  the  continuity  equation 
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(with  Z  measured  normal  to  the  plane),  so  that  if  we  write  R,  uR  as  a  dimension  and 
velocity  of  order  unity.  0(1).  then  Z,  are  each  0(5)  where  5  is  the  thickness  of 
the  shear  layer.  Since  terms  involving  squares  and  products  cf  fluctuating  quantities 
will  be  small  compared  with  those  involving  products  of  fluctuating  quantities  and 
mean  flew  derivatives,  the  only  terms  of  importance  in  (A-5)  are  for  R  »  0, 
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if  mean  velocity  derivatives  of  0(1)  are  neglected. 


(A-7) 


A-ii 


APPENDIX  B 


The  Evaluation  of  an  Integral 


It  was  shown  in  Section  2  that  the  pressure  correlation  depends  on  the  value  of 

-  _  3  /  1 

dr,  S„,r,  ^  — 


where  r  =  z  -  y. 

Kraichnan4  has  evaluated  a  similar  integral  following  the  method  of  Feynman23.  He 
found  that 


y  a.  p 

v 


=  2tt 


r  -l! 


where  -  x'  -  x  and  v,  the  volume  of  integration,  extends  to  infinity. 
If  in  (B-l)  we  put 
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then  with  ^(Sj,  0,  g3) 
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and  if  following  Feynman  we  introduce  a  new  variable 
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(B-5) 


and  following  Kraichnan4  we  use  the  identity 
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then 
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and 
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But  because  £7  =  0  we  see  that  s2  -  y2  and  so  from  (B-3),  (B-5)  and  (B-9) 
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Now  the  orders  of  integration  with  respect  to  r2  and  y2  can  be  reversed 
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equation  (B-10)  becomes 
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which,  apart  from  the  integration  over  r,  Is  just  half  the  value  found  in  (B-2)  above 
for  g  =  1,  when  the  integration  with  respect  to  y  is  taken  over  the  whole  Bpace  and 
not  the  half-space  considered  in  our  problem. 

The  error  introduced  by  the  use  of  (B-12),  in  place  of  the  exact  integral  (B-ll), 
is  known  to  be  small. 
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APPENDIX  C 


The  Structure  of  the  Large  Eddies  near  the  Wall  in  a  Turbulent  Boundary  Layer 


The  work  of  Townsend0 ■ 24  and  of  Grant10  has  shown  that  a  relatively  simple  structure 
exists  for  the  large  eddies  in  a  turbulent  flow,  or  at  least  the  experimental  observa¬ 
tions  are  not  inconsistent  with  the  hypothesis  that  the  large  scale  notions  can  be 
adequately  described  in  terms  of  relatively  simple  eddies.  Grant’s  measurements  of 
the  Rij  (r)  velocity  correlations  at  different  heights  from  the  surface  in  a  turbulent 
boundary  layer  give  rise  to  the  suggestion  that  the  large  eddies  near  the  surface 
[i.e.  in  the  constant  stress  layer)  have  the  form  of  ‘two-dimensional’  Jets  of  fluid 
originating  near  the  viscous  layer,  and  may  arise  from  the  instability  of  that  layer, 
and  are  roughly  aligned  in  the  direction  of  mean  motion.  Similarly  in  the  ‘outer 
region’  the  large  scale  turbulent  motion  again  appears  to  be  dominated  by  the  presence 
of  mixing  jets  of  turbulent  fluid  which  originate  in  the  interior  of  the  flow  and 
penetrate  to  the  region  of  the  non-turbulent  flow  outside  the  boundary  layer. 


A  simple  extension  of  the  Townsend-Grant  description  of  the  large  eddy  structure 
for  such  motions  led  us  to  assume  that  the  eddies  randomly  distributed  over  all  space 
may  be  described  by 
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where  (u.,  u?,  u,)  are  the  components  of  velocity  in  the  directions  (x , ,  x2,  x3),  with 
x,  in  the  direction  of  the  freestream  and  x2  perpendicular  to  the  wall,  and 
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The  velocity  components,  given  by  Equation  C-l  to  C-3,  satisfy  the  equation  of 
continuity 
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with  (a,,  a„,  a,)  having  Quite  arbitrary  values. 

i  2  3 


Following  Townsend  we  next  find  the  two-point  velocity  correlations  are  given  by 
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dy  (1  -  y2)  {l  -  (y3  +  r3)2}  exp(-  yf/2)  exp(-  (y  +  r)2/2)  (C-6) 


R22(r) 


and 


R33(D 
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dy  yj^!  +  ? j) (1  -  y3){l  -  (y3  +  r3)2)exp(-  y2/2)  x 

x  exp(-  (y  +  r)2/2) 

dy  y +  rj)(y3  +  ?3>(1  -  y 2> (1  -  y2  -  r2)  x 


(C-7) 


X  exp(-  y2/2)  exp(-  (y  +  r)2/2) 
where  y2  =  a1y1;  y2  =  a2y2;  y3  =  a3y3  and  dy  =  dy3  dy2  dy3. 


(C-8) 


After  integration  these  relations  can  be  written  in  the  following  form,  if  Rlx  etc. 
are  the  velocity  correlation  coefficients,  and  the  subscript  (b)  denotes  the  big 
eddy  contribution: 
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we  see  that 


(C-13) 


where  u,  etc.,  represent  the  big  eddy  contribution  to  u,  etc.  In  the  constant  stress 
b  * 

layer  Laufer  finds  that 


—r-  0.5  and 

u: 


whilst  in  the  outer  region 


=  =  0.5  and  ^  =  0. 2 


This  suggests  that  in  the  constant  stress  region 


=  2.5  a. 


a3  =  01  ai 


whilst  in  the  outer  region 
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=  0.5  a, 


These  results  show  clearly  that  the  role  of  the  u2  velocity  component  is  changed 
in  passing  from  the  constant  stress  region  to  the  outer  region. 

The  sijple  eddy  structure,  in  the  constant  stress  region,  portrayed  by  these 
results  is  that  of  an  elongated  vortex  ring  having  its  longest  direction  in  line  with 
that  of  the  mainstream.  It  is  rotating  in  a  plane  parallel  to  the  wall,  whilst  the 
longitudinal  motion  Is  wavelike,  being  away  from  and  then  towards  the  wall.  The 
latter  motion  is  not  unlike  the  jet-like  motion  described  by  Townsend  and  thus  may 
similarly  be  interpreted.  In  the  outer  region  the  osculations  are  diminished  in  scale 
whilst  the  eddy  spreads  out  slightly  in  a  direction  parallel  to  the  wall  and  perpend¬ 
icular  to  the  freestream. 

The  complete  value  for  R22  must  be  the  sum  of  the  contributions  from  the  small  and 
big  eddies.  Thus  to  R2J(r)b  we  must  add  on  (1  -  a2)  tines  the  value  for  R22(r) 
obtained  from  the  smaller  eddies.  If  we  assume  that  these  eddies  are  isotropic,  in 
contrast  with  the  strong  anisotropy  of  the  big  eddies,  we  have  that 


R22(r)  s 
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where  f(r)  is  the  longitudinal  velocity  correlation  coefficient  and  f'  5  df/dr.  Thus 
finally 


On  the  assumption  that 


R22  R22(s)  H  R22(b) 


f(r)  =  exp(-  r/I) 


(C-15) 
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the  values  of  R?7  have  been  evaluated  and  are  compared  with  Grant’s  results  in 
Figures  11,  12  and  13  for  the  measurements  made  close  to  the  wall  (i.e.  y/8  %  0.08). 
Although  the  agreement  is  not  good  quantitatively  we  have  qualitatively  obtained 
results  which  predict  the  basic  form  of  the  measured  results.  Similar  comparisons  can 
be  made  for  R13  and  R33  both  in  the  constant  stress  layer  and  in  the  outer  region 
provided  the  appropriate  values  of  ar  a2  and  a3  are  chosen  as  explained  above. 

With  these  values  for  R22(r)  we  can  now  find  the  two-point  pressure  covariance  on 
the  wall,  with  zero  time  delay.  Prom  the  previous  analysis  we  find  that  the  pressure 
covariance  on  the  wall,  with  the  separation  vector  is  given  by 
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and  the  corresponding  spectrum  function  is 
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Thus  if  P(x  ;  k\  fl)  is  the  scalar  wave  number  spectrum  function  on  the  surface  found 
by  averaging  tt  over  all  angles  in  the  plane,  noting  that  k  -  /k2^  +  k2  and 


P(x  ;  0  ;  0)  =  P(7  \  k  ;  0)  dx 


then 
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since  in  Equation  (C-18) 
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We  can  now  find  the  separate  contributions  to  P  from  the  small  and  big  eddies. 

Thus  with  R22(b)  Kiven  by  Equation  (C-10)  we  find  after  some  tedious  though  straight¬ 
forward  algebra  that 
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Similarly  with  R22(s)  given  by  Equation  (C-14)  with  f(r)  =  exp(-r2/i2)  (in  order  to 
simplify  the  algebra  but  not  necessarily  changing  the  order  of  magnitude  of  the 
results)  we  find  that 
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A  comparison  between  Equations  (C-22)  and  (C-23)  shows  that  whereas  the  small  eddy 
contribution  to  P  is  dominant  near  wave  numbers  of  k  -  l  (the  typical  energy  bearing 
scale  of  the  smaller  class  of  eddies)  the  big  eddy  contribution  is  very  dependent  on 
the  relative  scales  of  the  big  eddies  in  the  directions  (x^  x?,  x3)  respectively, 
i.e.  l/ar  l/a2,  l/a3. 

Now  in  obtaining  all  these  results  it  has  been  assumed  that  P-22(r)  does  not  change 
significantly  with  distance  from  the  wall,  which  is  certainly  not  true  if,  as  we  have 
shown  above,  we  pass  from  the  constant  stress  layer  to  the  outer  region.  A  rough 
comparison  of  Equation  (C-22)  with  Willmarth’ s  results6  clearly  shows  that  since  high 
frequency  peaks  beyond  «8  =  l (where  Sj is  the  boundary  layer  displacement  thickness) 

are  not  present,  the  pressure  spectrum  must  be  determined  from  the  form  of  the  velocity 
correlation  function  which  is  not  too  near  the  wall.  This  means  in  fact  that  a3/a 
cannot  be  large,  or,  in  other  words,  the  pressure  spectrum  is  largely  determined,  up 
to  the  high  frequency  cut-off  point,  from  the  big  eddy  contributions  with  eddies  whose 
transverse  extent  (in  planes  parallel  with  that  of  the  wall)  is  not  greatly  different 
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from  that  in  the  streamwise  direction.  This  would  imply  that  the  large  eddies  outside 
the  constant  stress  region  play  a  more  significant  role  in  the  determination  of  the 
wall  pressure  fluctuations,  than  do  the  large  eddies  inside  that  region.  This  is  not 
a  result  we  might  have  expected  from  considerations  of  the  region  just  outside  the 
viscous  layer  in  which  the  maximum  values  of  the  turbulent  intensity  and  mean  shear  act. 
However,  it  must  be  borne  in  mind  that  our  analysis  is  necessarily  very  approximate 
and  the  reason  for  this  result  might  appear  more  obvious  if  we  had  worked  through  the 
analysis  retaining  the  variable  mean  shear  and  turbulent  intensity  instead  of  replacing 
them  by  an  averaged  value  taken  over  the  entire  shear  layer. 

Finally  in  this  Appendix  it  is  worth  pointing  out  the  differences  which  exist 
between  our  results  and  those  of  Kraichnan5.  Kraichnan  makes  the  assumption  that  the 
turbulence  near  the  wall  is  homogeneous  in  planes  parallel  to  the  wall  but  not  in 
planes  normal  to  the  wall.  He  next  assumes  a  model  for  the  turbulence  in  which  ‘mirror¬ 
like’  velocity  boundary  conditions  are  satisfied  on  both  sides  of  the  wall.  In  this 
model  the  ut  and  u3  velucity  components  are  finite  on  the  wall  but  u2  is  zero.  The 
effective  correlation  coefficient  is  then  taken  as 

R22(x2,  x2;  rv  r 3)  =  R22(x2  -  x';  rr  r3)  -  R22(x'  +  x2;  r,,  r3)  (C-24) 

and  when  inserted  into  Equation  (C-17)  it  leads  to  the  result,  if  f  =  exp(-  a2  r2), 


Numerically  this  result  is  very  little  different  from  that  given  by  Equation  (C-23) 
but  is  considered  less  satisfactory  in  view  of  the  velocity  boundary  conditions  being 
essentially  different  from  those  existing  in  the  region  outside  the  viscous  layer. 
Kraichnan  also  includes  the  case  of  variable  mean  shear,  following  the  law 
t  ~  exp(-/3x2),  and  although  this  treatment  is  preferable  to  the  averaged  mean  shear 
approach  which  we  have  adopted,  the  results  are  not  qualitatively  different  and  do  not 
lead  to  an  order  of  magnitude  difference  in  the  numerical  results. 
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APPENDIX  D 


On  Space-Time  Correlations  of  the  Fluctuating  Velocity  and  Pressure 


In  the  analysis  above  the  pressure  covariance  has  in  the  main  been  calculated  for 
the  case  of  zero  time  delay  in  axes  fixed  in  the  wall.  However,  in  the  detailed 
analyaifl  the  structure  of  the  turbulence  (i.e.  R22(r))  is  required  in  a  frame  moving 
with  the  turbulence  and  it  is  not  immediately  obvious  what  the  relations  are  between 
the  turbulence  in  the  moving  and  fixed  frames.  This  problem  will  now  be  analysed 
below  ar.d  the  results  compared  with  the  experimental  results  of  Favre20  and  VfiTtmart.h6. 

Let  us  coiisider  a  field  of  homogeneous  turbulence  in  a  frame  of  reference  moving 
with  the  moan  velocity  up  past  the  fixed  point  P.  In  this  frame  of  reference  the 
four-dimensional  two  point  velocity  correlation  coefficient  is 

<■%  (r'> 

where  the  four  dimensional  vector  r ’  between  the  two  points  in  the  moving  frame  is 

given  by 


~r:  =  (r r '  r ' , 

1*  ?’  3’  '  u' 


and  r'  is  the  co-ordinate  representing  the  time  delay.  For  example,  if  the  longitudinal 
velocity  correlation  coefficient  were  given  by 

f(r')  =  exp(-  a2r'2) 


then  alternatively 


aV2  =  a2r'2  t  a2r/2  +  a2r'2  +  a2r'2 

1  1  ^2  2  “3  3  4*  4 


f(r')  -  exp(-  cr2r'2)  exp(-  v2r2) 


wIhmi1  1/a  is  a  measure  of  the  eddy  length  and  l/v  is  a  measure  of  its  life  time,  and 

r'2  _  r'2  +  r'2  +  r'2 

1  1  1  1  2  1  3 

This  result,  although  not  general,  is  typical  of  the  connection  between  space  and 
time  correlations. 

In  a  stationary  frame  of  reference  ( r )  in  which  the  turbulence  sweeps  past  at  the 
hpeed  u'jj  in  the  direction  rr  the  relation  between  r '  and  r  is 

with  r'  s  (rj,  r',  r';  t)  and  r  =  (r,,  r?,  r3)  , 


(r,  -J  up  (r')dr',  r2.  r^ 


We  now  make  the  hypothesis  that  the  turbulence  structure  is  unchanged  in  passing 
from  fixed  to  moving  co-ordinates.  Thus  with  up  independent  of  time,  the  velocity 
correlation  coefficient  in  fixed  areas 

(f)Rlj{r1  -  up  r,  r2,  r3;  r)  -  (n,)RiJ(f,>  (D-3) 

Examples 

(i)  Zero  time  delay  r  =  0 

(f)Rij(r1,  r2,  r3;  0)  =  (m)Ri:)(?';  0)  (D-4) 

This  means  that  the  spatial  correlation  of  velocities  in  a  fixed  frame  of  reference 
equals  those  in  a  moving  frame  of  reference. 

(ii)  Zero  spatial  separation  in  fixed  axes  r  =  0 

(f)Rjj  (-  up  r,  0,  0;  r)  =  =  -  up  r,  0.  0;  r)  (D-5) 

This  means  that  the  autocorrelation  in  a  stationary  frame  is  equivalent  to  a  space 
(r[  -  -  up  t)  and  time  (r)  correlation  in  moving  axes.  This  is  Taylor’s  hypothesis, 
which  is  supported  by  Pavres*  20  measurements  in  both  free  and  shear  flow  turbulence. 

(iii)  Optimum  time  delay  in  fixed  axes  Tj  =  Tm  up 

<f>R1J  (0.  0.  0;  rm)  (m)R1J  (0.  0,  0;  r  =  (D-6) 

This  shows  that  the  autocorrelation  using  an  optimum  time  delay  in  fixed  axes  is 
equivalent  to  finding  the  autocorrelation  in  a  moving  frame  of  reference. 

These  three  examples  are  made  clearer  by  reference  to  a  particular  velocity 
correlation  coefficient. 

Suppose  that  in  the  moving  frame 

(n,)Rn(r{,  0,  0;  r)  =  exp(-  a2  rj2)  exp(-  v2  t2)  (D-7) 

then  in  case  (i) 

^f)R11(r1,  0,  0;  0)  =  exp(-cr2  r'2) 

-  exp(  -  cr2  r2)  (D-8) 

showing  clearly  that  no  change  of  scale  occurs  in  changing  from  fixed  to  moving  axes. 

In  case  (ii) 


(f)RJ1(-  up  t,  0,  0;  t)  =  exp(-  o'2  u2  t2)  exp(-  v2  r2) 


(D-9) 


D-ii 


But  in  the  moving  frame  v  is  of  the  order  of  ct  u;  (where  uf  is  the  turbulent  intensity) 

so  that  if  u'  «  u„ 

p 

up  r,  0,  0;  t)  ^  exp(-  a 2  u2  t2)  (D-10) 

which  is  equal  to  (D-8)  if 


up  T 
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This  result  is  the  more  usual  form  of  expressing  Taylor’s  hypothesis. 


In  case  (Iii) 


<f)Rn(0.  0.  0;  rm) 


,  1 2 
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This  shows  that  the  autocorrelation  in  the  moving  frame  can  be  obtained  from  the  auto¬ 
correlation  in  the  fixed  frame  (by  comparing  (D-10)  with  (D-12))  by  replacing  the  time 
delay  (t)  in  the  fixed  frame  by 


T 


(D-13) 


The  results  of  Favre20  for  both  grid  and  wall  turbulence  (Fig. 14),  confirm  this 
simple  result  provided  that  the  time  (t  )  does  not  exceed  roughly  the  times  over  which 
the  autocorrelation  in  fixed  axes  are  essentially  non-zero. 


Similar  results  might  be  expected  to  hold  for  the  pressure  covariance.  Thus  if  we 
take  the  pressure  autocorrelation  in  moving  axes  to  be 

(n°P(0,0,0;  r)  =  (l  -  2v2  r2)  exp(-  r2)  (D-14) 

then  we  might  expect  that  in  moving  axes 


(ra)P(rJ,  0,  0;  t)  =  (l  -  272  r'2)  (1-2  v2t 2)  exp(-  v2  r[2)  exp(-  v2r2)  (D- 15) 

Therefore  in  case  (i)  above 

(f)P(rL,  0,  0;  0)  =  (1  -  r2)  exp(-  cr2  r2)  (D-16) 


In  case  (ii) 


<f)P(-  up  r,  0,  0;  r)  ^  (1  -  2cr?u2)  exp(-  a2 u2  r2)  (D-17) 

This  result  cannot  strictly  apply  to  the  pressure  covariance  at  the  wall  for  there 
Up  =  0.  However,  if  we  assume  that  the  turbulence  is  roughly  convected  over  the 
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plftne  wall  at  some  mean  speed  uc*  then  the  autocorrelation  in  fixed  axes  is  found 
from  (D-17)  with  u_  replaced  by  u„. 


In  case  (iii) 
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If  we  now  consider  the  general  case  of  a  space-time  correlation  in  fixed  axes  where 
r  t  t%%  the  optimum  time  delay. 
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where 


a  =  |r1/uc  t  -  1| 


According  to  Willmarth6,  from  measurements  of  the  pressure  fluctuations  on  wall 
turbulence  we  find  that 


uc  a*  0.8  Uj  (D-20) 

where  Uj  is  the  freestream  velocity.  A  comparison  between  Willmarth' s  results  and 
those  evaluated  from 

(f)P(rj  -  uc  r,  0,  0;  r)  -  (1  -  cr|r.  -  uc  r|)  (l  -  v | t | )  x 


x  exp( -  cr  |  Tj  -  ue  r  |  -  v  |t|)  (D-21) 

are  shown  in  Figure  15.  The  agreement  is  good  for  small  time  delays  (i.e.  up  to 
0.5  x  10~3  secs.,  which  is  roughly  the  time  for  non-zero  autocorrelation)  but  is  poor 
at  large  time  delays. 

The  reason  for  Willmarth’ s  results  being  non-negative  for  large  time-delays  probably 
results  from  the  growth  of  the  eddies  as  they  move  downstream.  For  small  time-delays 
this  effect  is  small  but  as  the  time-delay  is  increased  the  eddy-life  appears  to  be 
extended  an  a  result  of  this  growing  scele  of  the  moving  turbulence. 


•In  Incompressible  flow  the  tine  variation  of  the  pressure  in  fixed  axes  is  the  same  as  the 
time  variation  in  the  velocity.  Thus  the  autocorrelation  of  the  pressure  in  fixed  axes  must 
correspond  to  the  autocorrelation  of  the  velocity  averaged  over  all  angles.  If  then  the 
turbulence  has  a  dominant  convection  speed,  this  speed  will  also  apply  to  the  pressure  auto¬ 
correlation. 
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APPENDIX  E 


The  Evaluation  of  the  Pressure  Spectrum 


It  has  been  shown  above  that  the  pressure  covariance  at  the  point  x,  with  separa¬ 
tion  distance  g  cn  the  wall  is  given  by 


00  CD  00 

P(x;  0.  £3;  r)  =  <  t22  u2  >  J  dr2  f  drx  f  dr3  it2,(x;  r;  r) 

0  -oo 


r2  +  (r3  • 


k-rp 

If  the  velocity  correlation  coefficient  is  given  by 


(E-l) 


R22(r';  r)  =  e 


-v2t2  -  r ' 2/ i 2 


1  - 


r ' 2  +  r,'2 


1‘ 


(E-2) 


where  rr  is  the  separation  vector  in  co-ordinates  moving  with  the  turbulence  and  we 
assume  that  the  turbulence  is  convected.with  the  constant  speed  uc,  then  in  fixed 
co-ordinates 


R22(r;  t) 


2  2 
-V  T 


-(r2  +  rl)/l 2  -(r  -  u  r)2/!2 

e  2  3  e  1  c  x 


x  [l  -  r23/l2-(r1  -  ucr)2/l2]  (E-3) 


If  we  write  rt  =  r^l  etc  and  v  -  vl/ uc;  r  =  ucr/i  then 


P(?;  0;  r)  =  —  i2  <  r2  u2  >  e"v?  T? 
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u; 

5 

-J 
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(E-4) 


and  the  temporal  spectrum  function  P(x;  0;  <A  is  given  by 
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[l  -  r2  -  (F,  -  r)2]  dr 


where  o>  =  wZ/uc. 


On  taking  the  real  part  of  E-5  we  find  that  it  reduces  to 
P(x;  0;  w)  |  _  [  f  _  hi  +  r2 


4  P2  Is  <  t22  u|  >/n2uc 


-^2  +  r2  +  r2>) 

dr,  I  dr?  I  dr3  |  — 2— ^  ]  e\ 1  2  3/  x 
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x  e‘*v  +  1)t  cos(St)  [(l  -  r2  -  r2  -  T2)cosh(2  F,  f) 


+  2  rt  t  sinh  (2  r,  r)]  dr 


(E-6) 


Now  v  %  u'/uc  and  is  a  small  quantity  (i.e.  v  «  1).  We  can  therefore  neglect  the 
term  in  v2  in  Equation  E-6. 

Evaluation  of 

CO 

J  e'T  cos  &  r  [(l  -  F2  -  F2  -  r2)  cosh(2  F ,t)  +  2  F,r  sinh(2  F^)]  df  =  I, 
o 

This  integral  is  reduced  to  3  standard  integrals  which  can  be  evaluated  by  the  use 
of  Fourier  Transforms.  Thus 
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cos  cot  e~T  cosh(2  r  ,r)  df  (l  -  F2  -  r2) 
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Evaluation  of 
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From  Equation  E-7  we  find  that 
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Evaluation  of 
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If  we  put  X  =  Xt2  +  r2  and  we  use  polar  co-ordinates  (Y,  <p)  in  place  of  (r2,  r3) 
then,  from  Equation  E-8, 
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If  we  combine  the  above  results  then  from  Equation  (E-6) 
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We  note  finally  that 


2  u  - .  _ 
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(E-ll) 


which  is  exactly  double  the  value  found  by  Kraichnan"  for  isotropic  shear  flow  turbu¬ 
lence  without  boundaries,  and  provides  an  adequate  check  on  the  algebra  leading  to  the 
spectrum  function  given  by  Equation  (E-10).  To  allow  for  anisotropy  a  multiplying 
factor,  a  ^  1/3,  should  be  inserted  into  Equation  E-ll. 

The  spectrum  function  is  plotted  in  Figure  4,  where  the  abscissa  is  (oi^/IIj).  5^  is 
the  boundary  layer  displacement  thickness  and  Uj  is  the  freestream  velocity  outside  th® 
shear  layer.  If  we  take  uc  =  0.8  u,,  as  given  by  Willmarth6  and  Harrison7,  and 


E-i  v 
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I/Sj  =2.4  then 


coS  x 


o> 
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If  further  we  insert  into  (E-ll)  the  crude  estimate  (based  on  Laufer’ s  measurements 
of  t,2  and  y/ujand  evaluated  at  the  value  of  u7.t>/v  =  1550)  S^/u^,  <  r^2  >  =  2  x  10‘5, 
and  include  the  anisotropy  factor  of  a  =  1/3,  then 


P(s:  0) 


M  P 


2  4 


& 
Vi  P  K 


9. 3  x  10" 3 


compared  with  Willmarth’ s  value  of  0.006  and  Harrison's  value  of  0.0095  under  comparable 
conditions. 


Similarly  the  peak  value  of  P(x;  0;  co),  written  in  the  form 


2  P(x;  0;  co)  Ug> 
%  P2  ui 


10~u 


whereas  Willmarth  obtained  a  value  of  4  x  10"5  and  Harrison  a  value  of  6. 3  x  10"5. 


It  should  be  noted  that  our  estimates  are  only  3  db  higher  than  the  measured  values 
of  Willmarth  and  give  some  support  to  our  analysis.  The  inclusion  of  the  anisotropy 
factor  is  justified  by  Kraichnan  and  was  similarly  used  by  Lilley19  in  evaluating  the 
magnitude  of  the  pressure  fluctuations  in  the  mixing  region  of  a  jet.  Finally  it  might 
be  noted  that  better  agreement  between  theory  and  experiment  might  be  expected,  if  the 
anisotropic  form  of  the  velocity  correlation  function  were  used,  with  its  consequent 
improved  scale  relations  to  replace  the  ratio  1/8  as  used  above. 
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